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Brezin and Hikami ^21 El EJ El have considered a random Gaussian 
Hermitian ensemble with external source, 



1 -lTr(Af-A)^ M; 
Zj 

where M is random and A is deterministic. Notice this matrix ensemble, 
which had come up in the prior literature [2TJ EH EJ , ceases to be unitary 
invariant. The matrix A is chosen so that the support of the equilibrium 
measure has a gap, when the size of the random matrices tends to infinity. 
Through a fine tuning of A, the gap can be made to close at the origin. Brezin 
and Hikami propose a simple model having this feature, where the matrix A 
is diagonal, with two eigenvalues a and —a of equal multiplicity. Thus, upon 
letting the size of the matrix go to oo and after appropriate rescaling, they 
discover a new critical distribution, specified by a kernel involving Pearcey 
integrals [22J and having universality properties. 

P. Zinn- Justin [2U establishes the determinantal form of the correla- 
tion functions for the eigenvalues of the finite model. Then extending classi- 
cal connections between random matrix theory and non-intersecting random 
paths, Aptekarev, Bleher and Kuijlaars in [Zj give a non-intersecting Brow- 
nian motion interpretation of this Gaussian ensemble with external source. 
They also show that multiple orthogonal polynomials are the right tools for 
studying this model and its limit (see [HI El UH1 El) 

The present paper studies the Gaussian Hermitian random matrix ensem- 
ble 7i n with external source A, given by the diagonal matrix (set n = k\ + h-i) 



(a 



A :-- 



and density 



O 



\ 



o 



—a J 



_ e -1r$M*-AM) dMt 



I h 



(0.1) 



(0.2) 



Given a disjoint union of intervals E := Ui=i [&2i~i , &2i] C M, define the 
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algebra of differential operators 

& = (0-3) 
Consider the following probability: 

P n (a; E) : = P( all eigenvalues g£) = | / e - Tl ^ M ' 2 - AM) dM (0.4) 



In P], we have shown that for A = 0, the probability for this Gaussian 
Hermitian ensemble (GUE) satisfies a fourth- order PDE, with quadratic non- 
linearity (reducing to Painleve IV in the case of one boundary point) : 

(B\ + {An + 6B\ logP^B^ + 3B 2 - 4B_iBi + 6£ ) logP n = 0. 



The first question in this paper: Does the integral \0.4\ ), with A as in \0.1\) . 

satisfy a PDE? Indeed, we prove: 

Theorem 0.1 The log of the probability P n (a; E) satisfies a fourth-order 
PDE in a and in the endpoints b\, b<i r of the set E, with quartic non-linearity: 



(f+B^G- + F'B^gA (f+B^F- - F-B^F A 
- (F+G- + F~gA [f + B\F- - F-B\F + ^j = 0, (0.5) 



where 1 



F + := 2B_ 1 (^--B_ 1 )togP n -4A:i, 
oa 



2G+ ■= { H t>F + } B _ x -{H},F+} G- = G 



a — > —a 
k% «-» ^2 

a — > —a 
fci <->• ^2 



1 in terms of the Wronskians {/. g}x = gXf — JXg. 
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with 



H+ := |- (Bo - a|- - a£_i J logP n + f i5 i3_i + 4^- ) logP n 



&2 

+ 4A;i a + " 



ff+ := ^ (# - - aB_i^ logP n - (B - 2ai3_ x - 2) B_! logP n . 

(0.6) 

Remark: : The change of variables a i— > —a, &i «-> &2 hi the definition of _F~ 
and G _ act at the level of the operators. In fact, later, it will be clear that 
P n (a; E) is invariant under that change of variables. 

Again here we provide a natural integrable deformation of ()0.4j) (section 
1). As is well known, the probability for A = relates to the standard Toda 
lattice and the one- component KP equation (see pQ), the spectrum of coupled 
random matrices to the 2- Toda lattice and the two- component KP (see |2|), 
whereas we show that the model ()0.4|) relates to the three- component KP 
equation (section 2). This deformation enjoys Virasoro constraints as well 
(section 3), which together with the bilinear relations arising from 3-KP leads 
to the PDE of Theorem 0.1 (section 4). 

The second question concerns the Pearcey distribution, which we now 
explain. Following [Zj, consider n = 2k non-intersecting Brownian motions 
on R (Dyson's Brownian motions), all starting at the origin, such that the k 
left paths end up at —a and the k right paths end up at +a at time t = 1. 
As observed in |7j, the Karlin-McGregor formula [20J enables one to express 
the transition probability in terms of the Gaussian Hermitian random matrix 
probability P(a; E) with external source, as in (J0.4)) . 

Pq a (all Xj (t) e E) = lim I 

all 7 < - J E „ 
»li • • • . "fc — > — a 
S k +1, ■ ■ • , *2fc — * a 

1 

— det(p(t; 7i, ^))i<ij<n det(p(l - t; x v , <5p))i<t\j'<n 

(0.7) 




where p(t, x, y) is the Brownian transition probability 

1 



P(t,x,y) :-- 



(0.8) 



Let now the number n = 2k of particles go to infinity, and let the points a and 
—a go to ±oo. This forces the left k particles to — oo at t = 1 and the right k 
particles to +oo at t = 1. Since the particles all leave from the origin at t = 0, 
it is natural to believe that for small times the equilibrium measure (mean 
density of particles) is supported by one interval, and for times close to 1, the 
equilibrium measure is supported by two intervals. With a precise scaling, 
t = 1/2 is critical in the sense that for t < 1/2, the equilibrium measure 
for the particles is supported by one, and for t > 1/2, it is supported by 
two intervals. The Pearcey process V(s) is now defined as the motion of an 
infinite number of non-intersecting Brownian paths, just after time t = 1/2, 
with the precise scaling (see jjj): 



n 



2k-^ 



4. 



XjZ. 



th4 -+tZ 2 

2 



for z^O. (0.9) 



Even though the pathwise interpretation of V(t) is unclear and still deserves 
investigation, it is natural to define the probability 



¥{V{t) n E = 0) := limPj 17 " 2 (all x d (~ + tz 2 ) £ zE; 1 



5: j < Tl 
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Brezin and Hikami |T21 El El for the Pearcey kernel and Tracy- Widom 
[2Hj for the extended kernels show that this limit exists and equals a Fredholm 
determinant: 

W(V(t) n E = 0) = det (I-K tXB ), 
where K t (x,y) is the Pearcey kernel, defined as follows: 

p(x)q"(y) -p'(x)q'(y) +p"(x)q(y) - tp(x)q(y) 



Kt(x,v) 



where (note uj - 
1 

27 
1 

2tH 



p(x) 

q(y) 



id 

4 



p(x + z)q(y 



x-y 

z)dz, 



(0.10) 



e u 4 2 +u ydu = Im 



x 



due' 



4 



-uuy\ 
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satisfy the differential equations 

p'" — tp' — xp = and q'" — tq' + yq = 0. 

The contour X is given by the ingoing rays from ±ooe i7r//4 to and the 
outgoing rays from to iooe"" 1 "/ 4 , i.e., X stands for the contour 





/ \ 

The second result of this paper 2 is to give a PDE for the Pearcey distribution 
below in terms of the parameter t appearing in the kernel (jO.lOj) . Since this 
Pearcey distribution with the parameter t can also be interpreted as the 
transition probability for the Pearcey process, we prove: 

Theorem 0.2 For compact E = UI=i[ x 2i-i, #2i] and Bj = Yl?i x { +1 ~i^i 

Q(t;x 1 ,...,x 2r ) = log P (V{t) n E = 0) = log det(I-K tXE ) (0.11) 

satisfies a J^th order and 3rd degree PDE, which can be written as a Wron- 
skian 3 : 



l d 3 Q 

The proof of this statement, based on taking a limit on the PDE of 
Theorem 0.1 will be given in section 5. 



1 An integrable deformation of Gaussian ran- 
dom ensemble with external source 

Consider an ensemble of n x n Hermitian matrices with an external source, 
given by a diagonal matrix 

A = diag(ai, . . . , a n ) 



2 Tracy and Widom show in [23] the existence of a large system of PDE's involving a 
large system of auxiliary variables for Q and also for the joint probabilities at different 
times. 

3 given that {f,g}x ■= Xf.g ~ f.Xg. 
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and a general potential V(z), with density 

P n (M e [M, M + dM]) = — e- Tr{v{M) - AM) dM. 

Z n 

For the disjoint union of intervals E = IJI=i [hi-i, bzi], the following proba- 
bility can be transformed by the Harrish-Chandra-Itzykson-Zuber formula, 
with D = diag(Ai, . . . , A n ), 



P n (spectrum M C E) = ±- f e -K(V(M)-AM) dM 

Zk J Hn(E) 

Ai JE n 1 JU(n) 

1 r n 

= — A n (A)det[e-^^^] 1 < M < n r[rfA 4 



det[e a ^j] 
A n (A)A n (a) 



(1.1) 

In the following Proposition, we consider a general situation, of which (jl.lj) 
with A = diag(a, . . . , a, —a, . . . , —a) is a special case, by setting ip + = e az 
and tp~ = e~ az . Consider the Vandermonde determinant of the variables 
xi, . . .,x kl ,yi, ■ ■ .,Vk 2 , namely 

A n (x,y) : = A n (xi, . . -,x kl ,y h . . .,y k2 )- (1.2) 

Then we have 

Proposition 1.1 Given an arbitrary potential V (z) and arbitrary functions 
(p + (z) and tp~(z), define (n — k\ + k 2 ) 

(pi,...,p„) := e- y W (ip + {z),ztp + {z),...,z k ^ l ip + {z), 

l f-(z),ztp-(z),...,z k2 - 1 ip-(z)) . 

we have 
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„ n 

- / A n (z)det(pi(^))i<ij< n Y[dzi 

■ J E n ! 



ki 



k 2 



2- JE n 



A n (x,y)A fcl (x)A fc2 (y)J]^ + (^)e- yfe) ^n^~(^)e- V ' {K) ^ 



det 



z i+j - 1 < f + (z)e- v ^ 



E 



1 < i < fci 

< j < fci + k 2 - 1 



z l+j - l V -(z)e- y{ - z) 



Proof: On the one hand, using 

det (a ik ) 1<i k<n det (b ik ) k<n 



(1.3) 



1 < i < k 2 I 

< j < k t + k 2 - 1 / 



Yl det ( ai >°(j) h h°{j)) 



l<i,j<n 



a-eSn 



and distributing the integration over the different columns, one computes 



/n 
A n (z) det(pi(^))i< i j<„ JJ dzi 
i 



detO* ^i^ijxndet 



/ (4-V + (^)e- Vfe) ) i< 4 < fcl \ 

i < j < fci + k 2 



\ 



1 < j < fci + k 2 / 



= n! det 



^ +J '-V + (^)e" y(2) ^ 



E 



\ 



1 < i < fci 
0<j<fci+fc 2 -l 



^ +J '-V~(^)e" y(2) ^ 



ii 



1 < i < k 2 I 
0<3<fci+fc 2 -l / 



(1.4) 



On the other hand, one computes 

„ n 

j A n (z) det(pi(^))i<ij<„ TT ^ 



jEn *es n i=i 1 

„ n n 

^(-1) CT / A n (z a -i £),..., z a -i {n) ) JJpi(^i) Y[dzi 

jEn i=i 1 

„ n 

V(-1) CT / (-l) CT A„(^,...,z n ) JJpi(«j)dzi 

i=l 

fcl fc2 

A B (s, y) J] ar;-V + (*i)e~ V(iB0 ** II ^" VCtt)^^ 



A:i!A; 2 ! 



i 

fci fc 2 



1 1 



where A n (x,y) is defined in (|1.2|) . In the last identity, one uses twice the 
following general identity for a skew-symmetric function F(x%, . . . , Xk) and a 
general measure /i(dx), 



/ . . . ,Xk)A k (x) T\ v(dxi 

= / F(x 1 ,...,x h )j2(-ivn<tii^ dx J 

Jk " £s h 1=1 1 

= / 5Z( _i ) CTF ( xi '---' xfc )n^«M^ CT «) 

■ /R * ff GS fc i=l 

= / 5Z(-l) CT F(x (7 -i (1) ,...,x (T ^i (fc) ) JJx-~V(<fo.) 

J * h aeS k i=l 

= k\ F(xx,---,x k )T\x t i ~ 1 fi(dx i ). 



This ends the proof of Proposition 1.1. 



Add extra variables in the exponentials, one set for each Vandermonde 
determinant: 

t = (*i,t 2 , • • .)» s = (si,s 2 ,...), M = («i,«2, • • •) and /3 
Then, setting (n = fci + fc 2 ), 

2 00 

1 

Proposition 1.1 implies 

r klk2 (t,s,u;(3;E) := det m klM (t, s, u; (3; E) 

fcl &2 

fcl 



, „ fcl fc2 

4n / A n (x, y )n e sr^jj e Er^ 



Afci (a:) II e-^ +a ^' + ^e- -'Ur, 
j'=i 

A fe2 (y)n e 4-^-^ 2 e-Sr^ dy . ) j (1.5) 
j'=i 



where 



(nfj(t,s;P,E)) i<i< kl 
m fcl , fc2 (M,«;/3;E) := I { ^ U .^ E)) ^ 



< j < fci + fc 2 - 1 



< fc 2 

< j < fei + fe 2 - 1 



with 



H±(t,s;/3,E) = [ z ^^ e -^ +az+l3z2 e^^~ s ^ k dz 

J E 

Hj{t,u;P,E) = f z i+j - 1 e-£- az - f> * 2 eZ?(tk-uk)z k dz ^ (L6) 
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In particular, by (|1.3|) . the integral in has the following determinan- 
tal representation in terms of moments: 



and so 



?>,! 



A n (A) det[e" 



-V(Xj)+<HXj 



l<i,j<n 



( 



det 



z i +i- 1 e - 2 T- az dz 



\ 



1 < i < ki 

0<j<k 1 +k 2 -l 



1 < i < k 2 I 
0<j<k 1 +k 2 -l / 



P n (spec M C E) 



t=s=u=l3=0 



(1.7) 



(1. 



T klk2 (t,s,u;(3;M 

Remark: The integral enjoys the obvious duality: 

x < — > y, k\ < — > &2, t < — > t, s < — > u, a < — > —a, (3 < — > —f3. (1-9) 



2 Integrable deformations and 3-component 
KP 

Theorem 2.1 (Adler-van Moerbeke [5 ) Given the functions r nii „ 2 as in 
M.5)) . the wave matrix 

,(3)± \ 



W±, n2 {\;t,s,u) :-- 
with functions 

(A; *,*,«) : 



m(l)± ^(2)± 
v ni=Fl,n 2 V ™1>"2 



ni ,n 2 ±l 



(3)± 

ni=pl,n2±l 



(3)± 



ni,n 2 q=l ^ni±l,n 2 q=l *"lii2 



7"ni,n.2 (A ^0 



7"m ,n 2 \tj S, It) 
^T"2 c ±^M,A' 7 ""i."2(^ s ^' lt= F [A 1 ]) 
7~ni ,n 2 (£• S, If) 



(2.1) 
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satisfies the bilinear identity 



Wl k2 (\-,t,s,u)W e - ue2 (\-,t',s',u') T d\ = 
for all integers hi, ii, £2 > and t, s, u, t', s', v! G C°°. 



Proof: The moments, as denned in satisfy 



l'ij _ 



dh 



(2.2) 



9/4j 
ds k 

ds k 

and in matrix notation 



0. 



g/4 

du k 
duk 



dm r 



dt k 
where 



A := 



m A Tk 

//t OO,OO iS - 5 



/ 1 

10 

1 





dsi 



-A k m r 



(2.3) 



dui 



-A+moo^, (2.4) 



A O 
O O 



A_ 



O O 
O A 



Thus the moment matrix satisfies 



m r 



which implies the bilinear identity (J2.2j) . The details of proof can be found 
in |5] ■ ■ 
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Corollary 2.2 Given the above T-junctions T klk2 (t, s,u), they satisfy the bi- 
linear identities 



= i d\ \^-^e^-^r klM (HA- 1 ] ,s,u)r ei/2 (f +[A~ 1 ] , s', u>) 

J oo 

+ jd\ \^- 2 e^- s '^\ kl+1M [t,s-[\- 1 ] ,tt)r^_ 1A {t',s'+[\- 1 ] ,u') 
+ <fd\ \^~ 2 e^-<)y TkiM+l ( t , s , u -[A- 1 ] ) A _i (f, 5', u' +[A- X ] ) . 

(2.5) 

Upon specializing, these identities imply PDE's expressed in terms of Hirota's 
symbol 4 , for j = 1, 2, . . .: 

d 2 

S j(d t )r kl+1M or fcl _ ljfc2 = r fcifc 2 Q Sl Qt. +1 lo S Tfc i fc 2 ( 2 - 6 ) 

S j( 9 S )%-l,fa° 7 'W2 = T kik 2 Q tids . +i io & T kik 2 , (2.7) 



yielding 



d 2 \ogT kuk2 7- fcl+ i )fc2 rjfc 1 _ 1: 



A: 2 



dtidsx t? 



(2.8) 



fci,fc 2 



Alog^ = 4 SllQgTfcl ' fc2 (2-9) 

<9 2 1 _ 



lo g ^™ = ™« — (2.10) 
9si r *i-i.fe» a&7 lo g^i. 



Proof: The bilinear identity ()2.2|) yields nine identities, which are all equiv- 
alent, upon relabeling indices, to the tau-function bilinear identity ()2.5jl . 

4 Given a polynomial p(ij., £2, •••)> define the customary Hirota symbol p(dt)f o g := 
P^-Wi' ~By2' "•)/(* + 2/)s(* ~ 2/) ■ For later use, the s/s are the elementary Schur poly- 
nomials e£™ tiZ * := Y,i>o s i(t) zl and set s e (d) := s*(^, . . .). 
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Introducing standard shifts in the residue formulae 

t i — > t — a t' i — > t + a 
s i — > s — b s' i — > s + b 
u i — > u — c u' i — > u + c, 

and using Taylor's theorem, identity (|2.5j) is equivalent to 

oo 

^s, 1+£2 _ fcl „ fc2+J _ 1 (-2a)s J (4)e ErK ^ +fefc ^ +Cfe ^V £l , 2 o r felfe2 

3=0 
oo 

+ ^s, 1 _, 1+1+ ^-26)s J (4)e Er(afc ^ +6fe ^ +Cfc ^ ) r, 1 _ 1/2 or fcl+1 , fe2 

3=0 
oo 

+ £ s fe2 ^ 2+ i +J -(-2c)s i (4)e ErK ^ +6 ^ +Cfe ^ ) T 4A - 1 o r felife2+1 = 0. 

3=0 

(2.11) 

Taylor expanding in a,b,c and setting in equation (|2.11|) all a^fe^Cj = 0, 
except Oj+i, and also setting l x — k% + 2, £ 2 = &2> equation (|2.11j) becomes 

a j+l (- 2 Sj(dt)T kl+ 2,k 2 ° T kl ,k 2 + dg ^ dt , - r fci+l,fc2 r fci+l,fc2^) + = °> 

and so the coefficient of a J+ i must vanish identically, yielding equation ()2.6|) . 
setting fci — > fcj — 1. Setting in equation (|2.11|) all a^fe^Cj = 0, except fy+i, 
and £1 = hi, £ 2 = k 2 , the coefficient of in equation (|2.11|) yields equation 
(|2.7j) . Specializing equation (j2.6j) to j = and 1 respectively yields (since 
Si(£) = t± implies Si(d t ) = also s = 1): 



d 2 \ogT kuk2 T fcl+ i ifc2 r fcl _i. 



k 2 



dt i ds i t Im 

and 



& 



2 



kik'2 



d \ ( d 



7^-Tfcl+l,fc 2 J r fei-l,fc 2 — r fci+l,fc 2 I ^~" 7 >:i-l,*2 



Upon dividing the second equation by the first, we find equation ()2.9|) and 
similarly equation ()2.10|) follows from equation (|2.7jl . ■ 
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3 Virasoro constraints for the integrable de- 
formations 

Given the Heisenberg and Virasoro operators, for m > — 1, k > 0: 



_d_ 

dt n 



+ (-m)t- m + kS ,, 



j(2) 



3 2 



dUdtj 



2 E 



5 



i>l i+j=— m 



+ + 



m + 1 . / 9 



fc(Jfe + l) 



-m)t- 



2 ' V<9t 
we now state: 

Theorem 3.1 The integral r klk2 (t, s,u; (3; E), as defined in M.5)) satisfies 



B m r klM = V k ^ k2 r klM for m > -1, 
where B m and Y m are differential operators: 

2r „ 2r 

® m = b ? +1 gT> f° rE = U hi-i, hi] c R 
1 1 1 

+ al^H + (1 - 2/3)J^ 2 , fcl (-s) 



(3.1) 



and 



+lS fe2 (-«) - «J^i )fe2 (-«) + (1 + 2/3)J^ 2 , fc2 (-u) J 



(i) 



We state the following lemmas: 

Lemma 3.2 (Adler-van Moerbeke j3]) Given 



p = e v with — 



V 



9_ = E£A^ 
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the integrand 



dl n {x) := A n (x) Y[ (e^ hx *p(x k )d. 



k=l 



IX k 



satisfies the following variational formula: 



-^dl n (xi i-> Xi + ef(xi)xf +l ) 



(3-2) 

The contribution coming from Yii dxj is given by 



^a e (£ + m + l)j% tn dl n . 



(3.3) 



Lemma 3.3 Setting 

dI n = A n (x,y)f[e^^f[e^^ 

3=1 3=1 

kl x 2 

3=1 

& k2 {y O e ~ ~* ~ aV] e~^ u ^ dy 3 

3=1 

The following variational formula holds for m > — 1: 



—dl ( %i ^ Xi + £X 
de n \yi^yi + ey 



m+l 

i 

m+l 



(3.4) 



e=0 



Proof: The variational formula (J3.4)) is an immediate consequence of applying 
the variational formula (|3.2|) separately to the three factors of dl n , and in 
addition applying formula ()3.3|) to the first factor, to account for the fact 
that rijLi dxj Yij=i d-Vj is missing from the first factor. ■ 

Proof of Theorem 3.1: Formula ()3.1|) follows immediately from formula (|3.4jl . 
by taking into account the variation of BE under the change of coordinates. ■ 
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Using the identity, valid when acting on r klk2 (t, s, u; ft; E), 

d d d 

dt n 



ds n du n ' 

one obtains by explicit computation for m > — 1, 

T( 2 ) U\ _ ( m , 1HW 



v m 



+ «Jffi-i >fcl (-«) + (i - 2/3)J^ 2 , fcl (-^) 

I +JfS fa (-«)-^ 1 ^(-«) + (l + 2 ) 9)J« 2ifa (-«) J 



<9 2 



5 2 



<9 2 



\dtidtj dsidsj duiduj 



Ed .3 d 
I — 1" ls i— 1" m 



i>l 



i+m 



+ (h + k 2 ) ( + (-m)t-rn) - h ( t^— + (-m)s- m ) - k 2 ( \- (-m)u. 



dtr 



ds r 



d 



du r , 



771 ( 771 ~\~ 1 ) 

+(kl + hk 2 + kl)S m0 + a(h - k 2 )S m+lfi H (-t- m + s_ m + U- m ) 







+ a 



d d 
+ 



dt m+2 \ ds m+ i du m+ i 
d d 



+ (m + l)(s_ m+ i - w_ m+ i)^ 



+2/5 



rrt+2 <9 s m+2 



The following identities, valid when acting on T^kiifi s , u ] P] E), will also 
be used: 



d 
dsi 



1 ( d I jA _a_ _ 1 (_d_ i jA 

2 ^ da) ds 2 2 ^<9t 2 ~^ 8/3 J 

_d_ = _l ( _d a\ _9_ _ 1 f_d d_\ 

dm 2 \dt L da J du 2 2 \dt 2 8/3 J 

Corollary 3.4 The tau-function r = r kl)k2 (t, s,u; (3; E) satisfies the follow- 
ing differential identities, with B m = Y^T •' 
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V 9^_i 9si_i dui-i 
+a(k 2 - ki)r + {kiSi + A; 2 mi - (&i + k 2 )t 1 )r 



1 f„ 9 \ / 9 9 , 

- r = — + r 



2 V 9a/ \9si 9a 



isi- h m 



9.9 .9 
a,, . x 1 



i>2 



2 4^ V 9t,_i 9sj_i 9Mj_! 



+ 2 ~~ + -((fci + fc 2 )ti - fcisi - k 2 u 1 )r 



9 \ Or 
-[B -a—jT = — - (hi + k\ + k x k 2 )r 



i>i 

1 / 9 9 \ 9r 1 



„^9r >r-^ /. 9 .9 9 \ 



T 



9r 1 / 9 9 9 \ 



(3.5) 

Corollary 3.5 On £/ie /ocns £ = {t = s = -u = 0,/3 = 0} ; the function 
f = logrfc 1 fc 2 (t, s,u; j3; E) satisfies the following differential identities: 

|£ = -B^f + a(h-k 2 ) 



7 = K^-^) / + ^ 2 -^ 



9/ 
9s 



, j; v -i3 + a— ) / + k{ + hk 2 + k 2 
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df_ 

ds 2 

, d 2 f 
' ' dtidsi 

, d*f 
dt\ds 2 

, d*f 
' dt 2 dsi 



' ' da dj3 



(3.6) 



d_ 

da 

d d 



B-i)f- 2h 



da dp 



df 



a— + — - B + 1 - 2-^ - 2a(ifei - k 2 



da 



— (Bo ~a— + oB_0/ - B_i(Bb-l)/ - 2a{k x -k 2 ) (3.7) 



Proof: Upon dividing equations ()3.5|) by r and restricting to the locus C, 
equations ()3.6|) follow immediately. The essence of deriving ()3.7|) is that the 
Virasoro operators V n and the boundary operators B m commute. To derive, 
say, the first equation in the list ()3.7|) , rewrite the two first equations of ()3.5j) 

as 

df_ 
df 



-B-if 



+ a{k 2 -k 1 )+L 1 {f) + 



B-i 



d_ 

da 



f = ^ + ~a(h - k 2 ) + L 2 (f) + £ 2 
ds\ 2 



where Li are linear operators vanishing on £ and the li are functions vanish- 
ing on £,. This yields: 



d d 

di^^da- 



f 



(-B-J) 
1^/. d 



d 



d 



dsi-t 



dui-i 



-B-if) 



_d_ 

ds\ 



s^f. d . d . 
-> J I lU— \- ls i~ 1- lu i 



i>2 



d 



dti^x dsi-i dm-i 
y + (fcisi + k 2 u\ - (k\ + k 2 )ti) 



f 



d 2 



dsidti 



f + h. 
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The other identities (|3.7|) can be obtained in a similar way. 



4 A PDE for the Gaussian ensemble with ex- 
ternal source 



Proof of Theorem 0.1: First observe that, with n = k\ + k 2 , 

'H n {E) T klk2 (t,s,u;f3;R) 
An explicit computation over the whole range yields: 



F n (a;E) = ^- [ 

Z n JH, 



=u=/3=0 



r klk2 (t,s,u;p; 



t=s=u=l3=0 



hL AJX ' y) ( AkAx) W 

A fc2 (y)n 



j-+axi 



e 2 



-ay 



i=l 



c klk2 a klk2 e (kl+k2)a2/2 . 



This is obtained from the representation (jl.5|) in terms of moments, which 
themselves are Gaussian integrals, as shown in Appendix 1. From this for- 
mula, it follows that 



logr fclfc2 (t, s,u;/3; 



h + k 2 
= a 

t=s=u=0=O 2 



a 2 + k\k 2 log a + C, 



where c klk2 and C klk2 are constants depending on k±, k% only. It follows that 

hk 2 log a - C fclfc2 (4.1) 



logP n (a; E) = logr fclfc2 (0, 0, 0; 0; E) - ^!-^a 2 



Thus we need to concentrate on r klk2 (t, s,u; j3; E), which, by Theorem 2.1, 
satisfies the bilinear identity ()2.2j) and thus the identities ()2.9j) and (|2.1()jl of 
Corollary 2.2: 



9 . T k +i k 

— log — 

Oil Tkx-lfa 
dsi Tfn-x,k2 



dt2ds 



-\0gT klk2 



fclogr fclA;2 



dt 



a 2 



61 



(4.2) 
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whereas the first two Virasoro equations (|3.6j) yield, specializing to the locus 
£={t = s= u = 0,j3 = 0} and the indices ki ± 1, k 2 , 

9 . T k 1+ l,fc 2 Tfci+l,fc 2 . 

— log = -B-i log h 2a 

Otl Tjfci-l,jfc 2 r fc 1 -l,fc 2 

— log _^!_» = £_! - — log _ a (4.3) 
osi r kl _ lM 2 V oaj r kl -i,k 2 

From these three equations, the expression log Tfcl+1,fc2 can be eliminated, by 

first subtracting the first equations in ()4.2j) and ()4.3|) and then the second 

equations in (14.21) and (14.31) . Subsequently one eliminates log Tfcl+1 - fc2 from 

7 "fe 1 -i,fc 2 

the equations thus obtained, yielding 



or equivalently 

(4.4) 

Using the Virasoro relations ()3.7j) . one obtains along the locus £ = {£ = s 
u = 0,p = 0}: 

d 2 

4- 



dt\ds\ 

d 2 <)- \ . 9 



2 — — 2- 



dt 2 ds 1 dt\dsi 



logr fclfc2 =: 


F + 


logr fclfe2 = 


Hi 


J logTfcfc = 


Hi 



dt 2 d Sl dt x ds 2 J & ^ 1 ~ L dl3 

d 2 d 2 



2 «^T Pog^b = # 2 + (4.5) 
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where, using the identity (|4.1|) . along the locus £={t = s= u = 0,j3 = 0}, 

4fci 



F+ := 2B_x(— -B_ 1 )logn lfc2 -4fc 1 =2B_i(— -B-OlogI 



da 



d 



Hi := — [Bo- a- aB-i)logT klk 



da 



da 
d 



d 

B B-i + 4— ) logr fclfc2 + 2a(h - fe) 



= — I Bq — a— aB-i lo 



d_ 

da 
d 



da 
d 



r ( BoB^+4— ) log 



4afci 



4fcifc 2 



B -a—- aB-i ) logr felfe2 + (2aB_i-B + 2) B-i logr felfc2 + 2a(fc 1 + /c 2 ) 



3 



^ ,S - a^- - aB-i ) log! 



(2aS_i-S + 2)S_ilogI 



confirming ()0.6|) . Notice that the expressions above do not contain partials 
in P, except for the /3-partial appearing in the second expression of ()4.5|) . 
Putting these expressions (j4.5|) into ()4.4|) yields 



#-i^logr fclA;2 



1 ' 2 



2 , 2 



(4.6) 



and by involution a i— > —a, /3 i— > — /3, fci 
9 



where 



H7.-F- 
1 ' 2 



2 > 2 



=: G" 
,H7 



(4.7) 



a — > —a 
fei «-» fc 2 k\ «-» /c 2 

Remember the change of variables a i— > —a, /3 i— ► — /3, fci < — > &2 acts on the 
operators, since T^fa is invariant under this change; see (jl.6j) . 
Equations (|4.6j) and (|4.7|) yield a linear system of equations in 



Bj l ° g Z klk2 and B\ dl ° gTk ^ 



dp 



from which 



d log r fcl 



k-2 



dp 

G-F+ + G+F- 



dp 



, 91ogr fclfc2 



9/3 



G~(B_iF+) + G+(i3_iF- 
-F-(B^F+) + F+{B_ 1 F- 
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Subtracting the second equation from £>_i of the first equation yields the 
following: 

(f + B_iG~ + F'B^gA [f + B^F- - F'B^F 
-(f + G~ + F~gA [f + B\F- - F-B\F + ^j = 0, 
establishing Theorem 0.1. 



5 A PDE for the Pearcey transition proba- 
bility 

From the Karlin-McGregor formula for non-intersecting Brownian motions 
Xj(t), we have: 



P all x t (t) G E, l<i<n 



$iven Xi(0) = ji 
nven xAl) = 5, 



f 1 n 

/ — det(p(t;7i,x j ))i< i j< n det(p(l - t; x v , 6j>))i<i>,j><n J| dxi 

JE n 1 



1 (y-xf 

e t 



for the Brownian motion kernel 

P(t,x,y) :-- 



Aptekarev, Bleher and Kuijlaars introduce in [7j a change of variables trans- 
forming the Brownian motion problem into the Gaussian random ensemble 
with external source. For E : = UI=i[^2i-i, &2i] 5 we have, using this change of 
variables, 

in equality =, 

P^ a (all Xj(t) e E) 

all 2^(0) = 

P | all Xj(t) E E k left paths end up at —a at time t = 1, 
k right paths end up at +a at time t = 1 
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lim 

all 7, ; -> 

81, ... ,8 k —> -a 
<Wi> • • • 5 ^2fc — a 



__H i __2ay i 



/■ 1 n 

/ — det(p(t; 7i, Zj))i<ij-<n det(p(l - t; av, <V))i<i'j'<n JJ 

^ A n (x, y) (^A k (x) JJ e"^) + ^rfx^ f A*(y) J] 
in A n (x',2/) ( A fc (x')n e ~ % 

8=1 
' k 



i=l 



Fn{axl ^t ]bl \]to^y--- Mr p(i~t) 

with P n of Theorem 0.1, using (|1.1JI and (|1.3j) . with k = k\ = k%. Setting 



2t 



and e h{t) :-- 



t(l-t) 



2r 



5 



2r 



dv 



(5.1) 
(5.2) 



we find 



P± a (t; 6 2r ) = F n (ae^; he^, b 2r e h ^) = P n ( M ; v u v 2r ) • 

v=be ,h(t) 

(5.3) 

From Theorem 0.1, it follows that F n (u;vi, . . . ,v 2r ) satisfies the non-linear 
equation (|0.5j) . with a and all b^s replaced by u and Vi respectively. In order 
to find the equation for Pq a (t; 61, ... , b 2r ), one needs to compute the partial 
derivatives in t{ and b( in terms of partials in u and Vi, appearing in equation 
(10. 5 j) and use the relationship (|5.3|) . To be precise, compute 



0_ 

Of 



{BoYiB-xf Po a with i+j+i < 4 and i, j,£ > 0, (5.4) 
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yielding a system of 34 linear equations in 34 unknowns 

— j (B Q y (B^) e F n with i + j+£< 4 and i,j,£>0, (5.5) 

which one solves. Notice, one always writes (Bo) 3 (B-i) e in that order, using 
the commutation relation [B-i,Bo] = B-\. For instance, 

(B^yFp = ^{B^) j F n , (B yF± a = {B ) j F n j = l,...,4. 
" 1 g'(t)u^ + h'(t)B )P n 



<9t V 9m 







</(t) u — + fc'(t)£ </(*)«— + h'(t)B P n 



<9t 2 V 9m / V du 

+ ( <A*)«j^ + h"(t)B ) P„ 



The partials ()5.5|) thus obtained are now being substituted into the 4th order 
equation 1)0.5)1 . with a and 6j replaced by w and Wj, and thus the £>., by £>.,, 
yielding a new 4th order equation involving the partials ()5.4j) . 

Let now the number of particles n go to infinity, together with the corre- 
sponding scaling (see ff\ 123]) 

2 1 1 

n = 2k = — , ±a = ±— , 6j = XjZ, t = - + sz 2 , for z — > 0. (5.6) 

It is convenient to replace the ± in (J5.6)) by the variable e, which one keeps 
in the computation as a variable. The scaling combined with the change of 
variables ()5.3j) leads to the following expressions u and Vi in terms of z: 

2k = — A 

z 4 



u = ae 9 ^ = a 



2t eV2 /! + 



sz 



2 



1 — t z 2 V \ — sz 2 



h . e m = bj—^— = Xi (5.7) 
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So, the question now is to estimate: 

[f+B-xG- + F-B-iG+}(f+B-iF- - F-8-1F+ 
'f+G~ + F"G+) (f+B^F- - F-B 2 _ X F 



u 

Vi I ^ X 



%/2 / l+s* 2 



(5.8) 



For this, we need to compute the expressions F ± ,B-iF ± ,B 2 _ 1 F ± ,G ± and 
B_\G ± appearing in (|5.8|) in terms of 



' z {s;xi, . . . ,x 2r ) 
■■= logP 2/ , 
= Q(s;xi, . . . ,X2r) + O(z), 



1 


+ sz 2 




— sz 2 



; xi 



, X2r 



zV2 



3 - s 2 -* 4 



with 



(5.9) 
(5.10) 



. . .,x 2r ) = logdet (J - K s x EC ) , 

as shown in [213 . Without taking a limit on Q z (s; xi, . . . , X2 r ) yet, one com- 
putes 



F e 



1 



B\F E 



I6z 3 
1 

~32z 4 



16z 2 



32z 3 



2 9Q Z es 2 9Q 



Six 



9Q 2 



9s 16z 



9s 



0(z) 



0(1) 



3e 
8^9 



1 



128z e 



4z 7 

(B_i^)(B^iQ«) + mB B\Q z ~ {B\Q Z + 64s) B 2 ^ 



-%W 2 _ X Q Z + 16 



9Q Z 
9s 



9 3 Q 2 



9s 3 



+ 0(i) 



71^ 



36 Bi x Q ' £S Kl 



32z 10 
1 



16z 8 



512z' 



-(B.i^Lj^Q,) - 32B o B_i0, + (# 2 10 + 64s) 6^ 



-32S 3 X Q Z - 16B_i 



9 3 Q 2 



9s 3 



0(- 



9Q 2 
9s 

(5.11) 
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The formulae needed to obtain the expansions above for G £ and B_\G e are 
given in Appendix 2. From the expressions above one readily deduces 



- F~B-iF + = 



V2 
64z" 



-32(B - 2a4- - l)B\Q z 



\ 



-0(- 



e ^ B - lQz +0{- 4 ) 



V2z 



( 2{B^){B\Q Z ) 



F + B\F- - F~Bt\F + = 



F+G- +F-G+ = - 



2 IP+ 



1 

16z w 
4z 7 



-32(6 - 2sjL _ 2)B\Q Z 
V +{Bl 1 d -§f)(Bl 1 Q z )-l^ J 



+ %)■ 



Using this expressions, one easily deduces for small z, 

(f+B^G- + F~B-\G + ^j (f+B-iF- - F~B-iF+ 
-(F+G- + F~G + ) (F+B^F- - F-B 2 _ X F 



V-j i — * x 



'IT- 



1 _„2 
2 



2z 17 



( {B^,^ + (Bo-2)B^Q z } B _ \ 



+ 0(-Tb) 



£ 1 

= zr= ( the same expression for Q(s; x\, . . . , X2r) ) + O(-rr), 

2z L ' ' z w 

using (j5.10|) in the last equality. Taking the limit when z — ► yields equation 
(|Q~5T) of Theorem 0.2. ■ 



6 Appendix 1 

Setting 

/Xi+i-i(±a) :=^(*,s,«;/3,R)U =u= ^o= / ^"V^dz, 
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one computes 5 : 
Lemma 6.1 



T kl k 2 (t,s,u;[3,R)\ t=s=u=l3=Q = det 



/ (fj, i+ j(a)) o < » < fci -i \ 

< j < n - 1 



i {^i+j{ — a )) 0<i<fc 2 -l / 
\ 0<j'<n-l/ 



= c fclfe2 a fel/C2 e 2 a 2 . 

Cfclfc2 = (-2)^(2.)^ nYnV 



Proof: By explicit integration, one computes 

/io(a) = V2n e~ and /ij(±a) = v^r ( ±— ] e~2~ . 

\ eta / 

Define the Hermite polynomials (except for a minor change of variables) 

Pl (a) :=e~* (£) e<# = (± + ^(a). 

The following holds: 

P2i(a) = even polynomial, p 2 i+i( a ) = odd polynomial of a, 
which is used in equality = below, and 

Pk +n (a) = + M a)pt 1] + P2(a) P ( r 2) + ■■■+ PnPk, 
5 Remembcr n = k\ + k 2 . 
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where p k n ^ := (^)™£>fc and where Pi(a) are polynomials in a, independent of 
k; this feature is used in equality = below. Then we compute: 



T klk2 (t,s,u;[3,R)\ t=s=u=/3=Q 
= {V2^) n e~ det 



/ (Pi+j) o < i < fci - 1 

< j < n - 1 



V 



< j < n - 1 



(V2^) n (-1) — 2 — e~ det 



/ (Pi+j) o < i < fci - l 

< j < n - 1 



\ 



< j < n - 



= (V2^) n (-1) — 5 — e~ det 



/ 



V 



(p (l) ) 

Y 3 J < i < fci - 1 



\ 



< j < n - 1 



-iy P f 



< i < fc 2 - 1 ; 
< j < n - 1 / 



= c fcl fe 2 e 2 det 



1 < j < n 



l o,^fc 2 . 

\ 1 < 7 < n 



na 

= c fclfc2 e 2 det 



/ (a^V l ) i < i < fcj 

1 < j < n 



fc^ + 1 < i < n 

< j < n 



(i)-i+fei 



(tG5„ l<i<fei fci+l<i<n 

= c fclfc2 e^^(-l)^?(^)-)(a^) fe2 ^ 

l<i<n 

— C k!k 2 e a > 



where the ctjj are coefficients, some of which vanish. Indeed, each of the 
blocks in the matrix above is upper-triangular. To evaluate c' k k , observe, 



upon completing the squares in the exponentials and setting Xj 
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a, yj t— > i/j + a in the integral, 



T klk2 (t,s,u;{3; 



t=s=u=(3=0 

/ fel 2 



-TT7 / A fel+fc2 (x,?/) f A fcl (x)JT. 



fc2 «? 



A, ; , (//)]!' 1 "" ^ 

v i=i 

This integral equals 

= e (fel+fc2)a2/2 ((-2a) fclfc2 c feli oCo jfe2 + lower order terms in a) 

/ fci-i k 2 -i 

= e (fcl+fc2)a2/2 (-2a) fclfc2 (27r)^ j! JJ j\ + lower order terms in a 



o o 



The result in the first part of this proof implies the absence of the lower 
terms and thus Lemma 6.1. ■ 



7 Appendix 2 

In order to co mpu te the asymptotics 1)5.11)1 for the expression G £ and B-\G £ 
as defined in (jO.fij) . one needs the following asymptotics: 



1 „ e 6e 4 £ ,s 1 1 
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V2 4z 



9 W 1 S 2 dQz 1 B &2Qz 1 ( 4S 2 + 0^1 



^- li/2 = i B -^ +0{ ^ 



1 9 « m 1 r 3 dQz 6 r 2 1 r 93 4R 2 \k n + or 1 1 
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